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We present results on strange and multi-strange particle production in Au+Au collisions at 
^/sNN — 62.4 GeV as measured with the STAR detector at RHIC. Mid-rapidity transverse mo- 
mentum spectra and integrated yields of Kg, A, 3, f2 and their anti-particles are presented for 
different centrality classes. The particle yields and ratios follow a smooth energy dependence. 
Chemical freeze-out parameters, temperature, baryon chemical potential and strangeness satura- 
tion factor obtained from the particle yields are presented. Intermediate transverse momentum (pt) 
phenomena are discussed based on the ratio of the measured baryon-to-meson spectra and nuclear 
modification factor. The centrality dependence of various measurements presented show a similar 
behavior as seen in Au+Au collisions at ^/snn = 200 GeV. 
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PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

Strange hadrons in the final products of high-energy 
nuclear collisions provide valuable insight into the prop- 
erties of the created system, since they are not present 
inside the nuclei of the incoming beams. The produc- 
tion rates and phase space distributions of the strange 
particles may reveal different characteristics of the cre- 
ated fireball and thus they have been studied at different 
accelerators and experiments (H-Q. 

The enhancement of strangeness production in rela- 
tivistic heavy ion collisions was first proposed as a sig- 
nature of the Quark Gluon Plasma (QGP) in the early 
1980s [f|. The observation of an enhancement of strange 
baryon production relative to p + p collisions in SPS 
data [6j-|8|, confirmed later by data from the Relativis- 
tic Heavy Ion Collider (RHIC) @, has brought a lot of 
excitement to these studies. 

While the production yields are mainly determined 
by the yields in the low transverse momentum region 
(pr < 2 GeV/c), the intermediate transverse momentum 
region (j>t between 2 and 5 GeV/c) might give insight 
into the hadron production mechanism which in turn re- 
veals important characteristics of the system. One of the 
most interesting results observed in this region is the ap- 
parent enhancement of baryons over mesons which was 
originally observed in the p/n ratio but was also con- 



firmed in the strangeness sector with the A/ Kg ratio 11 



Different theoretical scenarios to explain such behavior 
are discussed showing that detailed experimental data 
are needed to discriminate between these proposals. 

We present in this paper results on the production of 
Kg, A, S and ft obtained by the STAR experiment at 
RHIC in Au+Au collisions at an energy of -y/ijvjv = 
62.4 GeV. This energy lies in the range between SPS 
and the highest RHIC energy and provides intermedi- 
ate measurement points in the energy excitation func- 
tions. The systematic study with detailed measurement 
of the excitation function of various observables is very 
important for understanding the QCD phase diagram of 
nuclear matter and for the search for a critical point of 
the phase transition that would correspond to the end 
point of the first order phase boundary as predicted by 

qcd mm- 

The paper is organized as follows. Section |H] is de- 
voted to the STAR experimental apparatus used for the 
measurements and the event selection. The analysis tech- 
niques to reconstruct the decays of the strange particles 
are described in section Mil along with the corrections 
needed to extract the physical spectra, the functional 
forms used to fit and integrate the pt spectra, and the 
procedure used to calculate the systematic uncertainties. 
In section lTVl the main results are presented and the phys- 
ical implications of these measurements are discussed, in- 
cluding detailed comparisons to the corresponding data 



from Au+Au at 200 GeV collisions |H|l|. Chemi- 
cal freeze-out properties are extracted from the particle 
ratios, strangeness enhancement is tested using the ratio 
of strange particles and pions, and baryon-to-meson dif- 
ferences are studied by comparing the Kg spectra to the 
A spectra. Final remarks and an overview of the general 
behavior of strange particle production considering the 
collision energy dependence are presented in section [V] 



II. EXPERIMENTAL APPARATUS 

The data presented in this paper were measured during 
the 2004 run of RHIC situated at Brookhaven National 
Laboratory (BNL) using the STAR detector. The data 
correspond to Au+Au collisions at ^/snn = 62.4 GeV. 
Descriptions and references concerning the STAR de- 
tector and its sub-detectors may be found in Ref. fl6j ] . 
The main subsystem of the STAR detector used for this 
analysis is a cylindrical time projection chamber (TPC) 
(length 4.2 m, inner radius 0.5 m, outer radius 2 m) cov- 
ering a pseudo-rapidity range of |?7| < 1.8 with a 2tt cov- 
erage in azimuth. The TPC is immersed in a magnetic 
field of 0.5 T parallel to the beam direction that enables 
the measurement of trajectories and momenta of charged 
particles with a transverse momentum (pt) above 0.15 
GeV/c. The TPC is filled with P10 gas (10% methane, 
90% argon) regulated at 2 mbar above atmospheric pres- 
sure. 

In addition to its track detection and momentum 
determination capabilities, the TPC provides particle- 
identification for charged particles by measuring their 
ionization energy-loss (dE/dx) in the TPC gas. Details 
and performance of the energy-loss particle-identification 
method are explained in Ref. [l7j ]. 

The collisions were triggered using a combination of 
several trigger detectors [181 ] . The minimum bias trigger 
requires coincidences between two zero degree calorime- 
ters (ZDC) located 18 m along the beam axis on each 
side of the interaction region. Beam-fragment (spec- 
tator) neutrons are not deflected by the beam guiding 
magnets and are thus detected in the ZDC detectors. 
In addition, a scintillation counter central trigger bar- 
rel (CTB) that surrounds the outer cylinder of the TPC 
was used to trigger on charged particle multiplicity in the 
mid-pseudo-rapidity region. Two beam-beam counters 
(BBC) were also used in the trigger system to measure 
the charged particle multiplicity in the forward region 
(3.3 < M < 5.0). 

The acquired data consist of a total of 6 x 10 6 minimum 
bias triggered events. In order to study the centrality 
dependence of strangeness production, these events were 
divided into centrality classes according to fractions of 
the total inelastic collision cross-section. The centrality 
of an event was defined as in previous analyses (l4l |. by 
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the measured (uncorrected) multiplicity in the TPC of 
charged tracks with specific quality cuts. 

In this analysis, only events with a primary vertex posi- 
tion within ±30 cm from the center of the TPC along the 
beam line were used to ensure a good rapidity coverage 
and relative uniformity of the detector efficiency. All the 
results quoted in this paper are mid-rapidity (\y\ < 1.0) 
results. For each event centrality class, the equivalent 
collision impact parameter is calculated based on a com- 
parison between the measured uncorrected charged parti- 
cle multiplicity dN ch /dr] and the calculated cross section 
using a Monte Carlo Glauber model [lj| . With this com- 
parison, it is possible to determine for each centrality 
class the average number of participant nucleons (N part ) 
and the average number of nuclcon-nuclcon binary col- 
lisions (-/Vbm)- These numbers together with the mid- 
rapidity charged particle yields dN c h/dy are shown for 
all centrality classes in Table HI 

III. ANALYSIS 

The analysis methods for the different particle species 
are detailed in the following sub-sections. The K s and A 
reconstruction method is described in IIII Al The recon- 
struction of multi-strange particles through the correla- 
tion of 3 particles is described in section UlIBI 

A. Kg and A reconstruction 

Kg, A and A were identified through the reconstruction 
of the weak decay topology: 

K s ->• 7T+ + n~ 
(branching ratio 69.20%) 

A — > p + 7T~ 

(branching ratio 63.9%). 

In addition to the dE / dx information of the daughter 
tracks, geometrical cuts such as the distance of closest 
approach (DCA) between the two daughter tracks and 
the pointing vector of the daughter tracks (distance of 
closest approach of the projected helix trajectory) away 
from the primary vertex position were used to improve 
the signal to noise ratio of the invariant mass peaks. A 
geometrical cut on the extrapolation of the reconstructed 
A trajectory from the primary vertex position was also 
used to reduce the contribution of feed-down A particles 
such as those from the 3 sequential decay. Cut selection 
was optimized considering a compromise between back- 
ground minimization and signal efficiency in the high-py 
bins. Table [TT| shows the geometrical cuts used in this 
analysis. 

The invariant mass distributions of the reconstructed 
particles, obtained after the geometrical and particle 



identification cuts, are presented in Fig. Q] panel (a) for 
Kg and panel (b) for A and A. The remaining back- 
ground underneath the mass peak was subtracted by 
making an interpolation of the spectrum on either side of 
the mass peak. Various polynomial functions were fitted 
to the background to estimate a systematic uncertainty 
for the signal area. The difference in the background es- 
timated using the two different methods was used in the 
calculation of the systematic uncertainty for each point 
in the final pr spectra. In general, the shape of the back- 
ground in the invariant mass spectra was smooth and 
relatively flat which resulted in uncertainties below 3%. 
The signal itself was then determined by counting the en- 
tries in the bins contributing to the peak and subtracting 
the estimated background. 

B. 3 and £7 reconstruction 

The multi-strange hyperons were reconstructed via the 
topology of their weak decays: 

5~ ->• A + 7r~ 

5+ -> A + 7T+ 

(branching ratio 99.887%), 

rr -> a + k~ 
n+ -> a + k+ 

(branching ratio 67.8%). 

A candidates reconstructed as described in the previ- 
ous section are combined with single track assuming it to 
be the bachelor 7r ± for the 3 decays and AT ± for the f2 
decays. 

Despite the identification of the daughters by their 
dE/dx, geometric selections were also used to improve 
the purity of the samples due to the large combinato- 
rial background. These selections were done by cutting 
sequentially on the different variables aiming for a com- 
promise between background minimization and signal ef- 
ficiency in the high-py bins. Table IIIII shows the geo- 
metrical cuts used in this analysis. The invariant mass 
of the assumed parent particle obtained from the candi- 
dates passing the selection cuts was calculated and plot- 
ted in a histogram, as shown in panel (c) for the 3 and 
(d) for the ft of Fig. Q] 

Even after the cut optimization, some background still 
remained in the invariant mass distributions. The back- 
ground under the invariant mass peak was estimated us- 
ing the same method as used for the A analysis, where 
polynomial functions were used to parameterize the back- 
ground outside the peak region. The background esti- 
mated by a polynomial fit function was subtracted from 
the invariant mass spectra in order to calculate the raw 
yield. In the following step, the raw spectra were cor- 
rected for the detector and analysis inefficiencies. 

C. Correction Factors 

These corrections include tracking efficiency and de- 
tector acceptance as well as reconstruction inefficiencies. 
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TABLE I: Summary of the centrality classes used to separate the events in this analysis. For each centrality class, the 
equivalent mean number of participant nucleons {N palt ) and the mean number of binary collisions are determined using the 
charged particle multiplicity distributions and a Monte Carlo Glauber model. Also shown in this table is the corrected charged 
particle multiplicity density (dN^/dy) that was obtained by adding the yields of protons, anti-protons, charged pions and 
charged kaons. Details of the analysis can be found in Ref. [141 ] . 





0-5% 


5 - 10% 


10 - 20% 


20 - 30% 


30 - 40% 


40 - 60% 


60 - 80% 


(iVpart) 


346.5 ± 2.8 


293.9 ± 4.2 


229.8 ± 4.6 


164.1 ± 5.4 


114.3 ± 5.1 


62.0 ± 4.9 


21.4 ± 3.0 


(iVbin) 


891 ± 57 


710 ± 47 


511 ± 34 


325 ± 23 


199 ± 16 


88 ± 10 


22 ± 4 


(dN ch /dy) 


582 ± 38 


476 ± 30 


359 ± 24 


249 ± 16 


166 ± 11 


86 ± 7 


27 ± 2 




1.1 1.12 1.14 1.3 1.35 

Invariant Mass [GeV/c?] 



1.75 



FIG. 1: Invariant mass distributions for selected candidates for (a) K s , (b) A, A, (c) E~, H + and (d) Cl + (note the different 
scaling factors for A, H and SI) for rapidity interval \y\ < 1.0 and pr > 0.2 GeV/c. The limits represented by the vertical 
dotted lines are the ones used for signal counting. 



The correction factors were determined as a function of 
Pt and the efficiencies for the different particles were cal- 
culated in different rapidity intervals between —1 and +1 
to verify that these corrections were independent of ra- 
pidity. The corrections were estimated by using Monte 
Carlo generated particles which were propagated through 
a TPC detector response simulator and embedded into 
real events on the cluster level. These enriched events 
were then reconstructed using the usual reconstruction 
chain. The correction factors for the different particles 
obtained by this method are shown as a function of pt 
in Fig. H 

The measured A spectra were also corrected for feed- 
down of weak decays by subtracting the contributions 
from the charged and neutral S decays. As already men- 
tioned in section IIII A[ a geometrical cut on the dis- 
tance of closest approach between the extrapolated re- 
constructed track of the A and the primary vertex posi- 
tion was used to reduce this contribution. Nevertheless, 
some of the secondary A still satisfied this cut, especially 
in the high-p^ region. The relative number of A parti- 
cles from the S decays that were in the primary particle 



sample was estimated using embedded Monte Carlo E! to 
calculate the reconstruction efficiency of these secondary 
A particles which was then scaled by the corrected yield 
of the measured S particles. The relative contribution of 
the feed-down was calculated for each pt interval and in 
the high-pr region was around 12%. Neutral S was not 
measured by the STAR experiment. Thus, to calculate 
the feed-down from these particles, the yield of the neu- 
tral S was considered to be equal to the yield of the ^ 
particles. The feed-down correction of the A spectra from 
the ft decay was not considered since it is negligible due 
to the lower yield of the f2 particles. The contribution 
from the Eo decay to the A yield was not subtracted as 
well due to the lack of S° yields measurements. 



D. Extrapolation of the pr spectra 

The final spectra for each particle species were ob- 
tained by dividing the raw yield pr distribution by the 
correction factors presented in the last section. 

As the spectra were not measured at low px due to 
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TABLE II: Geometrical cuts used in the analysis of Kg and A(A) production. 



Cut 


K° 


A and A 


DC A of VO to primary vertex 


< 0.5 cm 


< 0.8 cm 


DCA of VO-daughtcrs to primary vertex 


> 0.8 cm 


> 1.0 cm 


DCA between VO-daughters 


< 0.8 cm 


< 0.8 cm 


Number of hits in the daughters trajectory 


> 15 


> 15 


Radial decay length 


> 4 cm 


> 5 cm 



TABLE III: Geometrical cuts used in the analysis of the production of S, and their anti-particles. 



Cut 



and 



Q and SI 



DCA of parent to primary vertex 
DCA of bachelor to primary vertex 
DCA of V0 to primary vertex 
DCA of positive VO-daughter to primary vertex 
DCA of negative VO-daughter to primary vertex 
DCA between V0 and bachelor 
DCA between VO-daughters 
Number of hits in the daughters trajectory 
Parent decay length 
V0 decay length 
Daughter V0 inv. mass 



< 0.6 cm 

> 1.6 ■ -^/DCA of parent to primary vertex + 0.1 cm 

> 1.4 ■ -^/DCA of parent to primary vertex + 0.1 cm 



< 0.6 cm 

> 1.0 cm 

> 0.5 cm 

> 0.7 cm 

> 3.0 cm 

< 0.8 cm 

< 0.7 cm 
> 15 

> 5 cm 

1.1175 + 0.0125 GeV/c s 



< 0.6 cm 

< 0.6 cm 
> 25 

> 3 and < 20 cm 
> 2 and < 30 cm and > —4 * ^/"E decay length + 23 cm 
1.1150 + 0.0070 GeV/c 2 
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FIG. 2: Correction factors (acceptance x efficiency) for the 
most central events (0-5% for K%, A and S; - 20% for fi) 
at mid-rapidity (\y\ < 1) as a function of pr for the different 
particle species as obtained via embedding. The branching 
ratio of the measured decay channel is not factored into this 
plot. 



the limited acceptance coverage and at high px due to 
the finite statistics, the spectra needed to be extrapo- 
lated in these two regions in order to determine the total 
production yield (dN/dy) as well as the mean transverse 
momentum ((pt))- This extrapolation was done with 



two different functional forms, an exponential function: 
d 2 N 



2irpTdpTdy 
and a Maxwcll-Boltzmann function: 
d 2 N 



2TrpTdpxdy 



(1) 



(2) 



The spectra were better described with a Maxwcll- 
Boltzmann function, resulting in fits with lower \ 2 val- 
ues than the exponential function fits. In addition, the 
Maxwcll-Boltzmann function was also used in the anal- 
ysis of the data from Au+Au collisions at 200 GeV and 
thus allows for a consistent comparison between the two 
data sets. 

The main contribution for the determination of the 
dN/dy comes from the Iow-jjt region of the spectra. 
Therefore, in order to have a better estimate of the spec- 
tra shape in the low-p^ region the fit was performed con- 
sidering only the region of pr smaller than 1.5 GeV/c. 
This extrapolation varied from 5 to 30% of the final total 
yield for most particle species. In the case of the 11 par- 
ticles, the extrapolated yield was on the order of 50%, 
thus, the systematic uncertainty due to the extrapola- 
tion is large. The final dN/dy was obtained by integrat- 
ing the data in the measured region and using the fitted 
function in the low-p^ region. The contribution of the 
high-pT region to the dN/ dy was verified and found to be 
negligible. The difference in the total integrated dN/dy 
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using the Maxwell-Boltzmann and the exponential func- 
tion was considered to be the systematic uncertainty due 
to the extrapolation. 

The (pt) is obtained by integrating the whole range of 
the pr spectra as follows: 



(pi 



(3) 



The same functional forms used to determine the total 
yield were used to extrapolate the yield of the spectra in 
the unmeasured low-px region to calculate the {px)- The 
final uncertainty of the {px) values was estimated consid- 
ering the variation of the point to point uncertainties of 
the measured particle spectra, and the uncertainty in the 
contribution from the extrapolated low-p^ region. 



E. Systematic uncertainties due to the correction 
factors 

Several sources of systematic uncertainties were eval- 
uated and identified in the analysis. Variations in the 
detector performance and trigger efficiency during data 
taking and the dependence on position of the primary 
vertex were checked and found to be negligible. The 
main sources of the systematic uncertainties were the 
uncertainties in the correction factors of the measured 
spectra. 

In order to estimate the uncertainty due to the correc- 
tions of detector acceptance and efficiency of the anal- 
ysis cuts, a careful comparison was performed between 
the distributions of the geometrical parameters obtained 
from the embedded simulated particles and the distribu- 
tions of the tracks from real data. A good agreement 
between the data and the simulated embedded distribu- 
tions was obtained. The systematic uncertainties were 
calculated assuming a finite resolution to determine the 
cut value applied to the embedded sample during the 
correction factor calculation. The values of the cut were 
varied randomly given a certain assumed resolution and 
new correction factors were calculated. Next, the particle 
yields were obtained from these new correction factors. 
The variation of the yields obtained from this procedure 
was considered as one of the sources of the systematic 
uncertainty. This uncertainty was less than 10% from all 
centrality bins. 



IV. RESULTS AND DISCUSSIONS 

The corrected px spectra for K s , A, S, fi and their 
anti-particles are presented in Fig. [3j For better visual- 
ization, the spectra were divided by factors of 10, from 
central to peripheral data. The normalization factors are 
indicated in the figure. 



The shape of the spectra in the low-p^ region is sen- 
sitive to the characteristics of the evolution of the fire- 
ball such as transverse flow and the conditions of the ki- 
netic freeze-out. With the increase in the momentum, the 
shape of the spectra is also affected by soft jets, whose 
fractional contribution to the shape of the spectra in- 
creases with increasing momentum. The variations in 
the shape of the px spectra for different particles and 
different centrality classes can be better visualized by ex- 
amining the (px) values. The (pt) calculated from the 
data, using a Maxwell-Boltzmann function to extrapo- 
late for the unmeasured \ow-px region, is presented in 
Table H*Vl The difference in the calculated {px) between 
using a Maxwell-Boltzmann or an exponential function is 
presented as the systematic uncertainty. Figure 0] shows 
the evolution of {px) as a function of dN^/dy (see table 
U) for different particle species. 

The (px) values for A and S + are equal within uncer- 
tainties to the (px) for A and S~, respectively, and thus 
they are not plotted in Fig. 0J The {px) values from the 
f2 spectra are not included in this comparison due to the 
large uncertainties. The (px) for ir , K ± and p shown 
in Fig. [4] were presented in Ref. [l4[ • Error bars of each 
point correspond to the statistical uncertainty while the 
systematic uncertainty from the extrapolation is repre- 
sented by the gray bands. 

The (px) values of Kg are in agreement with the values 
of K~ and the A values with those of p spectra. The (px) 
shows a trend of increase from peripheral to central col- 
lisions for all particles. This increase is more pronounced 
for the heavier particles (p) than for the lighter particles 
(it). These observations are consistent with a collective 
radial flow among formed hadrons, which increases with 
centrality. 

The integrated yields for each centrality selection for 
the Kg, A and A are presented in Table|V]while the yields 
for the S~, S + , f2~ and + are presented in Table WT\ 
The systematic uncertainties correspond to the uncer- 
tainties due to the correction factors and the uncertainty 
due to the extrapolation. 

STAR experiment has presented K ± yields and spec- 
tra for the same colliding system in Ref. [3]. For the 
most central events, the integrated yields of Kg is lower 
than the average between the integrated yields of K + and 
K~ . Figure O shows the comparison of the Kg and K^ 
d 2 N / dydpx spectra in both top 5% and 30-50% centrality 
bins. The data points are consistent within uncertainties 
for the px range of 0.6-0.8 GeV/c but the Kg data points 
arc significantly lower for the two lowest px bins. The dif- 
ference is no larger than two sigmas. The charged Kaon 
spectra is limited by the dE/dx identification method to 
a maximum px of 0.75 GeV/c. In addition, the charged 
Kaon yields presented in Ref. [l4| are measured in the 
rapidity interval of \y\ < 0.1 whereas the Kg yields pre- 
sented are for rapidity interval of \y\ < 1.0. We also note 
that the difference between the charged Kaon spectra and 
the Kg spectra decrease with decreasing centrality and 
in the peripheral bin of 30-40% the data points are in 
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FIG. 3: Efficiency corrected pr spectra for the different centrality bins and for the various particles. Note that 7 centrality bins 
have been used for the K s and the A while only 6 and 3 have been used for the H and Q, respectively. Errors are statistical 
only. The A spectra are corrected for the feed-down of the 3 decay. 



TABLE IV: Average transverse momenta, (jot), in GeV/c, for the strange hadrons from Au+Au collisions at 62.4 GeV. The 
first error is statistical, while the second is the systematic error arising from the extrapolation in the low pr region. 







A 


A 




H + 





-5% 


0.684 ± 0.030 ± 0.023 


0.901 ± 0.041 ± 0.026 


0.885 ± 0.043 ± 0.023 


1.097 ± 0.037 ± 0.043 


1.130 ± 0.048 ± 0.050 


5- 


-10% 


0.694 ± 0.027± 0.023 


0.899 ± 0.039 ± 0.026 


0.892 ± 0.036 ± 0.025 


1.136 ± 0.040 ± 0.056 


1.136 ± 0.048 ± 0.050 


10 


-20% 


0.681 ± 0.024 ± 0.023 


0.902 ± 0.028 ± 0.026 


0.887 ± 0.026 ± 0.025 


1.115 ± 0.028 ± 0.049 


1.130 ± 0.037 ± 0.052 


20 


-30% 


0.660 ± 0.034 ± 0.023 


0.886 ± 0.032 ± 0.024 


0.892 ± 0.025 ± 0.026 


1.068 ± 0.025 ± 0.042 


1.042 ± 0.031 ± 0.035 


30 


-40% 


0.658 ± 0.036 ± 0.022 


0.862 ± 0.035 ± 0.021 


0.849 ± 0.029 ± 0.020 


40 


-60% 


0.627 ± 0.038 ± 0.022 


0.823 ± 0.035 ± 0.015 


0.807 ± 0.032 ± 0.013 


0.975 ± 0.039 ± 0.029 


1.022 ± 0.047 ± 0.041 


60 


-80% 


0.611 ± 0.062 ± 0.022 


0.772 ± 0.058 ± 0.010 


0.766 ± 0.075 ± 0.010 


0.869 ± 0.074 ± 0.013 


0.90 ± 0.11 ± 0.02 



agreement within the experimental error bars for all px 
points in the overlap region. 

Figure |6] shows the dependence of the singly-strange 
and multi-strange particle yields on the collision energy, 
including the results from this analysis. The lowest en- 
ergy points are from Au+Au collisions measured by the 
E896 experiment at the AGS [l| and the points in the SPS 
energy region are from Pb+Pb collisions measured by the 
NA49 experiment [2(| and the NA57 experiment . The 



RHIC data points for 130 GeV and 200 GeV 0,[H| are 
also presented. The strange baryon and meson produc- 
tion exhibits only a small energy dependence from AGS 
through SPS and on up to the top RHIC energies. From 
the AGS energy to the top SPS energy region, the A yield 
seems to increase at first, reaching a maximum around 
y/JJvw ~ 10 GeV and then starts to decrease slightly with 
energy. At RHIC the yield of A seems to stay approxi- 
mately constant, independent of the collision energy. 
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TABLE V: Integrated yield, dN/dy, for K s , A and A measured in Au+Au collisions at 62.4 GeV using data and a Maxwell- 
Boltzmann function for the extrapolation to the unmeasured low pr region. Quoted uncertainties are the statistical errors and 
the systematic uncertainties. The A and A yields are corrected by subtracting the contribution of the feed-down from the H 
weak decays. 







A 


A 





-5% 


27.4 ± 0.6 ± 2.9 


15.7 ± 0.3 ± 2.3 


8.3 ± 0.2 ± 1.1 


5- 


-10% 


21.9 ± 0.5 ± 2.3 


12.2 ± 0.3 ± 1.9 


6.1 ± 0.1 ± 0.8 


10 


-20% 


17.1 ± 0.3 ± 1.7 


9.1 ± 0.2 ± 1.3 


4.7 ± 0.1 ± 0.6 


20 


-30% 


12.1 ± 0.3 ± 1.1 


6.2 ± 0.1 ± 0.8 


2.99 ± 0.05 ± 0.40 


30 


-40% 


8.1 ± 0.2 ± 0.7 


4.1 ± 0.1 ± 0.6 


2.25 ± 0.04 ± 0.30 


40 


-60% 


4.0 ± 0.1 ± 0.3 


2.01 ± 0.04 ± 0.26 


1.16 ± 0.02 ± 0.16 


60 


-80% 


1.13 ± 0.05 ± 0.09 


0.504 ± 0.017 ± 0.07 


0.343 ± 0.012 ± 0.036 



1.4 



_n7I" 

oK" 

1.2h A P 



o 
> 

CD 
O 



Q. 

V 



1 

0.8 
0.6 h 
0.4 
0.2 




T A 

+ 

1 r> 



S 9 



□ □ □ □ □ □ □ 



10 



10 2 

dN cf /dy 



10 J 



FIG. 4: Extrapolated average transverse momenta (pr) as a 
function of dN c h/dy for different particle species in Au+Au 
collisions at 62.4 GeV. Statistical uncertainties are repre- 
sented by the error bars at the points while the systematic 
uncertainties are represented by the gray bars. The n, charged 
K and p data were extracted from Ref . [3| . 



Meanwhile, the fi and the anti-baryon yields show a 
clear and smooth rise with energy. The overall trend 
of the strange and multi-strange anti-particle yields with 
energy indicates a transition from a transport dominated 
production mechanism at the lower energies (AGS, SPS) 
to a production dominated by quark pair creation at the 
top RHIC energy. Transport dominated means that the 
initial baryon number of the colliding nuclei is trans- 
ported to the mid-rapidity region. This transition oc- 
curs in a way that the mid-rapidity baryon-free regime is 
reached smoothly. 

Figure [7] shows the anti-baryon to baryon ratio of A, 3 
and Q baryons as a function of energy which can be seen 
to reach almost unity for all particle species at y / ijvw = 
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FIG. 5: Kg dN/dpr spectra compared to the charged Kaon 
spectra for the event centrality of 0-5% and 30-40%. The 
charged Kaons data points are for rapidity range of \y\ < 0.1 
and were extracted from Ref. [Til ]. 



200 GeV. The dashed lines in Fig. [7] show the results of a 
thermal model calculation [22[ that will be discussed in 
section IW~A"1 

Figure [8] shows the anti-baryon to baryon ratio of p, A 
and 3 baryons as a function of the centrality of the col- 
lisions, represented by the measured number of charged 
hadrons per unit of rapidity (dN^/dy). Results from 
■y/sjviv = 200 GeV collisions (open symbols) and from 
the y?SNN = 62.4 GeV collisions (solid symbols) are com- 
pared. At t^snn = 62.4 GeV, the p/p and A/A ratios 
show a similar trend, decreasing with increasing central- 
ity, consistent with a higher degree of baryon number 
transport from beam rapidity for central collisions. The 
decrease with centrality seems to be larger at 62.4 GeV 
than at 200 GeV, which is consistent with a higher frac- 
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TABLE VI: Integrated mid-rapidity yield, dN/dy, for H _ , H + , and Cl + measured in Au+Au collisions at 62.4 GeV (using 
a Maxwell-Boltzmann function for the extrapolation). Quoted uncertainties are the statistical errors and the systematic 
uncertainties. 









q- 


0+ 





-5% 


1.63 ± 0.09 ± 0.18 


1.03 ± 0.09 ± 0.11 






5 


-10% 


1.16 ± 0.06 ± 0.16 


0.86 ± 0.08 ± 0.12 


0.212 ± 0.028 ± 0.018 


0.167 ± 0.027 ± 0.015 


10 


-20% 


0.96 ± 0.04 ± 0.11 


0.59 ± 0.03 ± 0.06 






20 


-40% 


0.52 ± 0.02 ± 0.06 


0.357 ± 0.016 ± 0.037 


0.056 ± 0.006 ± 0.005 


0.038 ± 0.005 ± 0.003 


40 


-60% 


0.183 ± 0.008 ± 0.021 


0.116 ± 0.005 ± 0.017 


0.0165 ± 0.0027 ± 0.0014 


0.0103 ± 0.0020 ± 0.0010 


60 


-80% 


0.042 ± 0.003 ± 0.005 


0.032 ± 0.003 ± 0.004 
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FIG. 6: Strange particle production yields at mid-rapidity 
in central Au+Au and Pb+Pb collisions versus the center of 
mass energy ^/snn- The top panel shows results for Kg and 
A. The AGS values are from E896 [H (centrality 0-5 %). 
The SPS values are from NA49 (centrality - 7 %) and 
the RHIC values are from STAR [J, [lf| (centrality - 5 %). 
For the multi-strange baryons E and Q (bottom panel), the 
SPS results are from NA57 [l(centrality 0-11 %) and the 
RHIC values are from STAR Hlllll (centrality 0-20 %). 



FIG. 7: Anti-baryon to baryon yield ratios for strange baryons 
versus the center of mass energy ^/sjvjv. A/ A is shown in the 
top panel while the multi-strange baryons are on the bottom 
panel. The data from AGS are not corrected for the weak 
decay feed-down from the multistrange baryons while the data 
from SPS and RHIC are corrected. The lines are the results of 
a thermal model calculation (see text section llV A[l . The AGS 
values are from E896 [H (centrality - 5 %). The SPS values 
are from NA49 H (centrality - 7 %) and the RHIC values 
are from STAR [J, []]| (centrality 0-5 %). For the multi- 
strange baryons 5 and Q. (bottom panel), the SPS results are 
from NA57 [Ifcentrality - 11 %) and the RHIC values are 
from STAR [ll,[2H (centrality 0-20 %). 



tion of net-baryons being transported from beam rapid- 
ity to mid-rapidity at the lower collision energies. At 
sjvjv = 200 GeV where there is a wider rapidity gap 
between the beam particles, the net-baryon density at 
mid-rapidity is lower, thus, the anti-baryon to baryon ra- 
tios of the studied particles are higher than at 62.4 GeV. 
In addition, the ratio is independent of centrality at 
^Jsnn — 200 GeV, indicating that pair production is 
dominating and that there is a smaller variation of net- 
baryon density with centrality. 



A. Chemical properties 

The chemical properties of the bulk particle production 
can be addressed by using thermal-statistical models (23r - 
[27j . They are applicable under the assumption that the 
yields of particles are governed by statistical laws. The 
particle abundance of species i per system volume (Ni/V) 
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FIG. 8: Anti-baryon to baryon yield ratios for strange parti- 
cles and protons as a function of dN ch /dy at ^/sjvjv = 62.4 
and 200 GeV. The p data were extracted from Ref. 0. The 
\/snn = 200 GeV strange hadron data were extracted from 



can be parameterized by 
9i 



Ni 
V 

where 



(2tt) s 



7 S ' 'exp( — )±1 



-I -1 



T 



cli 



d 3 p (4) 



(5) 



and the gi factor is the spin-isospin degeneracy, T c h is the 
chemical freeze-out temperature that marks the end of 
inelastic interactions, Hq, Hb> l^s are the charge, baryon 
and strangeness chemical potentials respectively, Bi, Qi 
and Si are the baryon, charge and strangeness numbers, 
Ei is the energy of the particle, 7 S is the strange quark 
phase-space occupancy and is the number of valence 
strange and anti-strange quarks in the particle i. 

These thermal-statistical models have been successful 
in reproducing the particle ratios for a large range of en- 
ergies and also have been able to reproduce ratios includ- 
ing strange and multi-strange particles such as the Q/w 
ratio (28l . [29| . The evolution of the anti-baryon to baryon 
ratio of the studied particles with ^snn is also well re- 
produced in these models as indicated in Fig. [7] by the 
dashed lines. These lines correspond to the anti-baryon 
to baryon ratios using a thermal model calculation with 
parameters (temperature and chemical potentials) given 
by the systematic study presented in Ref. [3C| . The data 
at the intermediate 62.4 GeV energy are also well repro- 
duced. 

In order to perform a thermal fit to the STAR data, 
we used the model THERMUS described in [H|. The 
parameters considered in this thermal fit are the tem- 
perature (T c h), the baryon chemical potential (/is) and 



1.2r— 

0.8 - 
0.6- 



fitted 



CD 

DC 0.4- 



0.2- 



ji,K, p 

7t, K, p, A, S 

Data points 



> 

Q 



Q. 



< 
l< 



1w 



l< 



+ o in 

IM ^ 



FIG. 9: (color online) Particle yield ratios as obtained by 
measurements (black dots) for the most central (0 — 5%) 
Au+Au collisions at 62.4 GeV, and statistical model predic- 
tions (lines). The ratios indicated by the dashed lines (blue) 
were obtained by using only 7r, K and protons, while the ra- 
tios indicated by the full lines (green) were obtained by also 
using the hyperons in the fit. 



the strangeness saturation factor (7 S ). The strangeness 
chemical potential (fi s ) was constrained by the initial 
S/V ratio (strangeness per volume). The charge chemical 
potential (/xq) was fixed to zero and the system volume 
was not considered because particle ratios were used. To 
calculate the particle ratios, proton yields were corrected 
to subtract the contributions from the weak decays of A 
particles. The feed-down contribution from the £ decay 
into the proton and anti-proton yields were done within 
the thermal model and the correction was on the order of 
7%. The pions were corrected for the weak decay feed- 
down, muon contamination and background pions pro- 
duced in the detector materials. The A yields were cor- 
rected for the weak decay feed-down from the 2. Figure[9] 
shows the particle ratios as obtained by measurements for 
the most central data as black dots. 

The fit procedure was performed with two different 
data sets, considering first only the particle ratios of p, 
p, 7T^~ and K^. Results from this fit are shown in Fig.[S]as 
dashed lines (blue online) and it is seen that the model 
fits quite well the ratios of pions, protons and charged 
kaons, but it under predicts the ratio The ther- 

mal fit was repeated considering also the particle ratios 
with A, A and E ± . The f2 ratios were not included in 
these thermal fits because the yields were obtained for 
different event centrality bins due to limited statistics in 
data. The results of the thermal fit, including the A and 
5 ratio are also shown in Fig.|H]by the full lines (green on- 
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FIG. 10: Chemical freeze-out temperature T" c h (panel a) and 
strangeness saturation factor 7 S (panel b) as a function of the 
mean number of participants. 



line) . The standard deviations of the data to the thermal 
fit result are also shown at the bottom of Fig. |9] Most 
of the results are within one standard deviation from the 
data for the case where the strange particle ratios are 
also considered in the calculation. The inclusion of the 
strange baryons in the thermal studies results in a more 
constrained fit and a better description of these yields. 
In addition there is a small increase of T c h from (147 ± 
7) MeV to (163 ± 5) McV but the value of 7 S remains 
constant at 0.88 ± 0.06. The particle ratios including the 
strange baryons are well described by the thermal fit, but 
it over-predicts the production of p with respect to pions. 

Figure \TU\ shows the evolution of the chemical freeze- 
out temperature (top panel) and the strangeness phase- 
space occupancy factor (bottom) with the number of par- 
ticipants. The solid symbols represent the results for 
tJsnn = 62.4 GeV and the open symbols represent the 
result for ^Jsmn = 200 GeV. T c h is independent of the 
centrality of the collision and its value is consistent with 
the top RHIC energy measurements. In contrast, 7 S in- 
creases from a lower value for peripheral collisions to 
saturate at a value consistent with unity for the most 
central events at y/ls~NN = 200 GeV and close to unity 
(0.87 ± 0.07) at ^/sj^ = 62.4 GeV. These results are 
similar to the results obtained previously at the higher 
RHIC energies, where 7 S rises from peripheral collisions 
and reaches values consistent with unity for the most 
central events [32j . 

Figure [Til shows the energy dependence of the chemical 
freeze-out temperature and the baryon chemical poten- 
tial. The open symbols were extracted from [22l |. The 
solid stars show the values obtained from this analysis 



FIG. 11: Temperature and baryon chemical potential ob- 
tained from thermal model fits as a function of y^ijviv (see 
[22| ). The dashed lines correspond to the parameterizations 
given in [%| . The solid stars show the result for ^snn 
and 200 GeV. 



62.4 



for Ts^FaF = 62.4 and for the 200 GeV STAR data with 
the same particle ratios. These results are consistent with 
other thermal fits presented in Ref. [22j for the same en- 
ergy range, where data points from other RHIC experi- 
ments are considered. The dashed curves in Fig. [IT] are 
the numerical parameterizations of T c h and /1 b presented 
in Ref. [30[ using: 



1.308 GeV 



1 + 0.273 GeV- 1 ^/! 



(6) 



and 



T, 



ch 



0.166 GeV 



-0.139 GeV _2 /4- 



0.053 GcV" 4 Mb 



(7) 



The values obtained in this analysis for T c h and fiB 
from the Au+Au 62.4 GeV data are in good agreement 
with the parameterizations that were obtained fitting the 
data points at higher and lower energies. As the collision 
energy increases, the baryon chemical potential decreases 
continuously reaching very small values at RHIC. In con- 
trast, the chemical freeze-out temperature (T c h) seems 
to achieve saturation at approximately 160 MeV and all 
three RHIC energy points shown in Fig. [TT] are already 
at this temperature. 



B. Strangeness Enhancement 

Strangeness enhancement is traditionally defined as 
the ratio of the strange particle yield measured in heavy- 
ion collisions, normalized by the mean number of par- 
ticipant nucleons ((N part )) and the yield measured in 
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p + p collisions. Strangeness enhancement was observed 
at the lower SPS energies Q as well as at higher RHIC 
energies Q. The multi-strange baryons show a higher 
enhancement factor than the A particles, which is con- 
sistent with the picture of the enhancement of ss pair 
production in a dense partonic medium as opposed to a 
pure hadronic medium It is also argued that part 
of the enhancement can be due to canonical suppression 
effects in p + p collisions [H|]. In this scenario, the pro- 
duction of strange hadrons is actually suppressed in p+p 
collisions due to limitations of phase space. To discrimi- 
nate among the different mechanisms that can affect the 
strangeness enhancement factors, it is necessary to study 
the systematic behavior for different collision energies, 
systems and centralities. 
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FIG. 12: Ratio of baryon (solid symbols) and anti-baryon 
(open symbols) to tt + as a function of dN c h/dy for ^/snn = 
62.4 GeV (left panel) and y/J^F = 200 GeV (right panel). 
The 7r and p data were extracted from Ref . [l4| • 

Since there are no data on strange baryon production 
for p + p collisions at ^/snn = 62.4 GeV for comparison, 
we show the ratios of strange baryons and tt mesons in 
Au+Au collisions as a function of dN^/dy normalized to 
the ratio in the most peripheral centrality bin (Fig. IT2"j) . 
The data points were not included in this plot due to 
the difference in the centrality selections. The measure- 
ments at y/sNN = 62.4 GeV (left panels) and 200 GeV 
(right panels) are shown for strange baryon to pion ra- 
tios in the top panels, and for anti-particles in the lower 
panels. The pht and p/n ratios have been included for 
comparison (I4| . 

It is seen that the strange baryon to pion ratio in- 
creases with centrality and the relative enhancement is 
higher than observed in the p/n ratio. The multi-strange 
baryon S is the particle with the highest enhancement. 
In the 62.4 GeV data, there are indications of a difference 
in the enhancement rate between the strange baryons and 
anti-baryons, with the latter showing a smaller increase 



with centrality. This difference can be due to the non- 
zero net-baryon density which makes it easier to form 
strange baryons than strange anti-baryons. In addition, 
the lower net baryon density at 200 GeV reduces the 
difference between strange baryons and anti-baryons at 
this energy. At 200 GeV, there is no significant difference 
between baryon and anti-baryon enhancement for all par- 
ticle species, but the hierarchy of enhancement observed 
at 64.2 GeV is maintained at the top RHIC energy. From 
the traditional strangeness enhancement data, where the 
yields are normalized by the p + p data, the onset of the 
enhancement is pronounced and is already present in the 
most peripheral heavy-ion collisions, consistent with the 
effect that would be expected from the canonical suppres- 
sion hypothesis. In this context, the increase observed in 
Fig. 1121 which is normalized by the peripheral bin of the 
measured Au+Au collisions, does not include the effects 
of the canonical suppression and could be exclusively due 
to the enhancement of the ss pair production mechanism. 
However, the p/n ratio also increases with the system 
size, reflecting the increase of the baryon stopping as the 
collisions becomes more central. This increase of the net- 
baryon density would also favor the production of A and 
5 with increasing centrality. 

Our results show that there is a net increase of the 
strange baryon production with respect to the pion yields 
as the system goes from peripheral to central collisions, 
suggesting a true increase of the ss pair production rel- 
ative to the light quark production in central collisions. 
This conclusion is further strengthened by the observa- 
tion of the enhancement in the <f>/ir ratio measured in 
central Au+Au collisions at 62.4 GeV with respect to 
peripheral collisions as presented in Ref. [HI, HH . 

Figure [TBI shows the ratios of strange baryons to ir me- 
son yields as a function of collision energy. The anti- 
baryon to 7r meson ratios increase as a function of energy 
while the A/ir and S/-7T ratios show a decrease from SPS 
to RHIC energies and then remain constant in the RHIC 
energy region. The trend in the energy dependence of 
the A/ir ratio was already observed in the K + /tt + ra- 
tio where the peak structure at the SPS energies is more 
pronounced [14| . This shape is consistent with the net- 
baryon density dependence on the collision energy [i"4| . 
The A/-7T and S + /7r ratios show a strong monotonic in- 
crease with energy. The Sl~/7r and & + /ir ratios show a 
smaller increase with energy, however the experimental 
uncertainties are large. The dashed lines in this figure 
correspond to the ratios of these particles using a ther- 
mal model calculation with parameters (temperature and 
chemical potentials) given by the systematic study pre- 
sented in Ref. [30| . The general trend is well reproduced 
by this calculation. 



C. Nuclear Modification Factors 

One of the most important results from RHIC was the 
suppression of the particle spectra in the high-py region 
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FIG. 13: Ratio of baryon (solid symbols) and anti-baryon 
(open symbols) to n~ as a function of ^/snn- The lines are 
the results of the thermal model calculation (see text section 
|IVA|). The SPS values are from NA49 [U (centrality 0-7 
%) and the RHIC values are from STAR [J GH (centrality 
— 5 %). For the multi-strange baryons 3 and fl (bottom 
panel), the SPS results are from NA57 @] (centrality 0-11 
%) and the RHIC values are from STAR pj, [U (centrality 
0-20 %). 



observed via the nuclear modification factor known as 
Raa- It is the ratio between the measured particle spec- 
tra from Au+Au collisions normalized by the correspond- 
ing number of binary collisions and the measured spec- 
tra in elementary p + p collisions !36] . Together with the 
measurement of the disappearance of back-to-back high- 
Pt hadron correlations in central Au+Au collisions at- 
tributed to interaction of the energetic jet particles with 
the formed medium (37j . it was concluded that matter 
created at RHIC is highly interacting. The identified 
particle Raa from Au+Au collisions at 200 GeV shows 
that even the strange particle spectra are also suppressed, 
and there is a grouping of baryons and mesons with re- 
spect to their suppression in the intermcdiate-pr region. 
This suppression seems to persist at lower energies and 
it was observed for pion, proton and charged kaons in 
Au+Au collisions at 62.4 GeV It is important to 

study whether strange baryons follow the same trend as 
at 200 GeV and whether the suppression also shows a 
baryon and meson grouping (37j . 

For the strange baryons measured at 62.4 GeV it is 
not possible to evaluate Raa due to the lack of p + p 
data. However, a comparison between the spectra mea- 
sured in central collisions and that from the peripheral 
collisions can also be used to analyse the suppression at 
high-pT. This comparison is done with the differential 



FIG. 14: Nuclear modification factor Rcp, calculated as the 
ratio between 0-10% central spectra and 40-80% peripheral 
spectra, for it, K g , A and S particles in Au+Au collisions at 
62.4 GeV. The n Rcp values were extracted from Ref. [Tot ] . 
The gray band on the right side of the plot shows the uncer- 
tainties on the estimation of the number of binary collisions 
and the gray band on the lower left side indicates the uncer- 
tainties on the number of participants. 



nuclear modification factor Rcp which is the ratio be- 
tween the central and peripheral spectra scaled by the 
number of binary collisions ((Vbi n )): 



Rcp 



d 2 N 



central 



/dp T dy 



d 2 NP eii P heTBl /dpTdy 



V 



peripheral 



a /"central 
iV bin 



(8) 



The result of the Rcp measurement comparing the 
spectra of central 0-10% data with the peripheral 40-80% 
data for the different particle species is shown in Fig. [T4l 
For comparison, the Rcp of pions from Au+Au collisions 
at 62.4 GeV is also included 0. 

The Kg data, presented as solid circles show that the 
Rc p reaches a maximum value of approximately 80% at 
Pt ~ 2 GeV/c and decreases to a ratio of approximately 
40% for higher transverse momentum. The general shape 
of the Kg Rcp curve shows the same trend as the n Rcp 
(shown as open circles). In the intermediate-pT region, 
between 1.5 and 4 GcV/c, the K s ratio is higher than 
the 7T ratio. The 7r Rcp reaches a maximum around 
0.7 at a pr value of approximately 1.5 GeV/c. But, in 
the high-pT region, both K s and 7r seem to be equally 
suppressed, saturating at around 40%. The difference 
between the Kg and ir curves at intermediate-pT is also 
observed in Au+Au collisions at 200 GeV [H|. 

The strange baryons, A (open squares) and 5 (solid 
triangles), show a different behavior. The ratio increases 
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FIG. 15: Nuclear modification factor Rcp, calculated as the 
ratio between 0-5% central spectra and 40-60% peripheral 
spectra, for A and S particles measured in Au+Au collisions 
at 62.4 GeV. The gray band corresponds to the equivalent 
Rcp curve for the A particles measured in Au+Au collisions 
at 200 GeV fil. 



to a maximum above unity at approximately 2.5 GeV/c 
and then seems to decrease towards unity without show- 
ing significant suppression up to 4.5 GeV/c. This be- 
havior is different from that of the proton Rcp which 
reaches unity at around 2 GeV/c and starts to decrease 
for higher values of px [HI • 

In Fig. [TH the Rqp curves of A and S measured in 
Au+Au collisions at 62.4 GeV are compared to the Rcp 
of A from Au+Au collisions at 200 GeV presented in 
Ref. 
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To allow for a direct comparison to the 200 GeV 
data, the Rcp was calculated from the ratio between 
spectra of centrality bins 0-5% and 40-60%. The shape 
of the Rcp curve for 200 GeV is similar to the results 
from 62.4 GeV, also showing a flat region after reach- 
ing the maximum, and the onset of suppression only 
occurs at higher values of pr- The maximum in the 
intermediate-j»T region seems to be higher at 62.4 GeV 
than at 200 GeV. However, considering the systematic 
uncertainty of the normalization factor, shown as the 
gray band in the right side of Fig. \TE{ it is not possi- 
ble to conclude that the ratio measured at 62.4 GeV is 
higher than the ratio measured at 200 GeV. 

In summary, the K s spectra at 62.4 GeV are consistent 
with a suppression scenario for px > 4.0 GeV/c while 
the A and 5 baryons seem to show no suppression for 
Pt ~ 5 GeV/c. In Fig.[14]one observes that the difference 
between the A and S curves is smaller than the gap to 
the Kg data in the intermediate-pT region, suggesting a 
baryon-to-mcson separation. However, at 62.4 GeV the 



baryon meson separation is not as clear as observed in 
the results from Au+Au collisions at 200 GeV Hal. 



D. Baryon to Meson Ratio 

Another interesting result related to particle spectra 
and yields obtained from RHIC data is the behavior of 
the proton to pion ratio with px, where an anomalous 
increase of baryon over meson production was observed 
in the intermediate-pT region and furthermore this in- 
crease was higher in more central collisions [Is} . This 
result is referred to as the baryon anomaly, and it is the 
focus of theoretical discussions that include particle pro- 
duction and also particle interaction with the medium. A 
baryon-to-mcson enhancement was also observed in the 
strangeness sector, where the ratio A over Kg also showed 
an anomalous increase in the intermediate-pT region [Til ] . 

Figure [TC] shows the centrality dependence of the A to 
K s ratio as a function of px- The ratio increases as a 
function of px reaching a maximum around 2.5 GeV/c. 
Above this px there is a slope change in the curve, and 
the relative meson production starts to increase with re- 
spect to the baryons. In the intermediate-pr region, be- 
tween 1.5 GeV/c and 4 GeV/c, there is a dominance 
of baryon over meson production. Also, the maximum 
value of the ratios diminishes when going from central to 
more peripheral collisions. One of the proposed theoreti- 
cal explanations for this behavior involves the parton co- 
alescence and recombination mechanisms in hadron pro- 
duction [40l - [42l | , where a baryon at a certain px results 
from coalescing partons of a lower px- In this scenario, 
the baryon spectra would have a higher average px than 
mesons, resulting in a ratio that increases with px- In 
general, the coalescence models have been successful in 
describing the increase of the baryon-to-mcson ratio, in- 
cluding the ratio with strange particles. 

An alternative proposal comes from the consideration 
of QCD higher-twist contributions in the initial state of 
hadron production [43||, where baryons would be pro- 
duced in a color singlet configuration, and thus would 
be color transparent and not suffer final state interac- 
tions, increasing the final yield of baryons over mesons. 
Within this proposal, the increase of the baryon-to-meson 
ratio with the increase of centrality would be naturally 
explained by the higher degree of meson scattering as the 
volume of the medium increases. 

In the high-pT region, the baryon-to-meson ratio starts 
to decrease with the onset of the contribution from jet 
fragmentation, in which meson dominates over baryon 
production. 

In order to discriminate between different theoretical 
models it is necessary to study the systematic behavior 
for different energies and different system size. Figure IT71 
shows the maximum value of the measured A/ Kg ratio 
for the different centrality bins (quantified here by the 
average number of participants) in Au+Au collisions at 
62.4 GeV and 200 GeV [ll[ as well as p + p at 200 GeV 
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44j . It is clear that the baryon-to-meson ratio exhibits 
the same trend at 62.4 GeV as observed at 200 GeV: an 
increase of the maximum with the increase of the cen- 
trality and the same ratio of approximately 1.7 and 1.0 
for the most central and peripheral bins, respectively. 
Figure [17| also shows the A/K s ratio for the 62.4 GeV 
data. The anti-baryon to meson ratio peaks at much 
lower value, below unity, but shows the same trend of in- 
crease with centrality. The difference observed between 



the A and A enhancement over the K s can be attributed 
to the non-zero net-baryon density. 



V. CONCLUSIONS 

In this paper we present the strangeness production in 
Au+Au collisions at ^/snn = 62.4 GeV as measured by 
the STAR detector at RHIC. For each particle and anti- 
particle which includes K s , A, 5 and f2, the yield and px 
spectra at mid-rapidity have been extracted for different 
collision centrality classes from peripheral to the most 
central events. 

Analyzing the anti-baryon to baryon ratio of the stud- 
ied particles as a function of y/sJ^N demonstrates that the 
mid-rapidity baryon-free regime is approached smoothly 
from the AGS energy to the top RHIC energy. The data 
at y/SNN = 62.4 GeV presented here fit well into the sys- 
tematics of the excitation function of these ratios. Also 
the production rates of all three anti-baryons studied 
show a smooth increase with energy whereas the baryon 
yield excitation functions are susceptible to the interplay 
between the increase of the pair production mechanism at 
mid-rapidity and the decrease of the net-baryon density 
with increasing rapidity range in higer energy collisions. 

The chemical composition of the formed system was 
studied in the framework of a statistical thermal model, 
where the strangeness saturation factor parameterizes 
the degree of equilibration between the strange quarks 
and the lighter u and d quarks. Results using this model 
to fit our data show that even for the lower energy 
62.4 GeV data, the system created in central Au+Au col- 
lisions is still consistent with a saturation of strangeness 
production, similar to the result obtained at 200 GeV. 

Another similarity with the measurements at 200 GeV 
is seen in the behavior in the intermediate-pT region. Re- 
sults of the baryon-to-meson ratio and the nuclear modi- 
fication factor point to a different behavior of mesons and 
baryons in this px region. Various theoretical proposals 
were discussed which try to explain these experimental 
observations. To test any of these models, or even to 
place tighter constraints on the models, it is necessary 
to perform a careful and systematic comparison over a 
wide range of energies. Within this context, the results 
presented here are a valuable complement to the exist- 
ing systematics and fill the large gap between the lower 
energy SPS data and the top RHIC data. 

As an overall conclusion, it seems that strangeness pro- 
duction in Au+Au collision at ^snn = 62 A GeV is 
qualitatively very similar in all aspects to the produc- 
tion at y/SNN = 200 GeV. All excitation functions vary 
smoothly from AGS to top RHIC energies. 
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